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We address the important question of how to characterize the normal state of fermionic superfluids 
under the influence of a strong effective magnetic field, implemented through rapid rotation or novel 
artificial field techniques. We consider the effects of crossing from BCS to BEC and the role of non- 
condensed pairs, or pseudogap effects. Using a simple extension of Gor'kov theory we demonstrate 
how these pairs organize above the transition T c into precursors of a vortex configuration, which are 
associated with distortions of the ideal Abrikosov lattice. This non-uniform normal state appears 
to enable "Bose condensation" in a field which is otherwise problematic due to the effective one- 
dimensionality of Landau level dispersion. 



The role of non-condensed pairs in fermionic and 
bosonic superfluids under the influence of a strong ef- 
fective magnetic field has not been clearly addressed in 
the literature. One could imagine that in anticipation 
of the vortex configuration of the superfluid phase, these 
non-condensed pairs exhibit some degree of inhomogene- 
ity. The goal of the present paper is to investigate this 
normal phase. Our study should lead to a reconsideration 
of previous work on rotating cold Fermi gases [1 where 
the upper critical rotation frequency was computed un- 
der the assumption that non-condensed pairs were not 
present, even in the BEC regime. It may also be rel- 
evant to "normal state vortex" models [2j [3] which are 
argued to be important in high temperature supercon- 
ductivity. Furthermore, we note the community excite- 
ment over proposals to establish strong effective magnetic 
fields through "artificial" means. These are implemented 
through either asymmetric tunneling in a lattice [4] or 
geometric gauge potentials [5j [6] . This makes our work 
on the interplay of rapid rotation and superfluidity par- 
ticularly topical. 

Finally, we believe our work bears on a puzzling aspect 
of superconductivity in high magnetic fields, even within 
the BCS framework. As has been emphasized in the lit- 
erature, the degeneracy of Landau levels implies that the 
fluctuations around the BCS phase are effectively one- 
dimensional [7] , representing the free propagation of par- 
ticles along the field direction; this is well known [8] to 
be problematic for stable superconductivity. 

In this paper we suggest that the origin of this latter 
difficulty is due to the fact that in BCS theory, Cooper 
pairs only exist in a condensed state. Here we show 
that by introducing non-condensed pairs associated with 
a finite pairing gap (or "pseudogap") at the onset of 
condensation, a strict "dimensional reduction" [7 is no 
longer present; three-dimensional behavior leading to sta- 
ble condensation at T c can then occur. Our approach is 
to be contrasted with previous attempts [9] to address 
the early concerns raised by Schafroth [8 in Bose gases. 



Here we emphasize the contrast with BCS theory, where 
the only pairs under consideration are in the condensate. 

One goal of this paper is to stimulate experimental 
searches for the predicted precursors of the vortex con- 
figurations below T c . Here we characterize the degrees of 
freedom associated with non-condensed pairs and demon- 
strate that they create pair density inhomogeneities. Al- 
though these excited pairs become progressively more im- 
portant as the superconductor crosses from BCS to BEC 
(and pairing becomes more stable), they may well be 
necessary even in the BCS limit for stable condensation. 
These ideas relate to earlier work on pair density waves 
[10], but are in contrast to previous studies on the inter- 
play of a pseudogap and magnetic field [11 where density 
inhomogeneities were not contemplated. 

To understand these non-condensed pairs we are 
guided by Landau-Ginzburg (LG) theory where in zero 
magnetic field non-condensed pairs are associated with fi- 
nite center of mass momentum q; these represent gapless 
excitations as q — >• 0. In non-zero field, the natural coun- 
terpart should be associated with slightly distorted con- 
figurations of the vortices. These pair excitations become 
gapless as they approach the lowest energy superconduct- 
ing vortex configuration. However, we also emphasize 
that these non-condensed pairs are distinct from previ- 
ously considered Landau-Ginzburg vortex lattice fluctu- 
ations which only address the condensate [I2j [13]. We 
now proceed to describe how we incorporate these non- 
condensed pairs, and their effects at and above the su- 
perfluid transition. 

Rewriting the Gor'kov Equations We begin with 
a BCS-type theory, and proceed to separate out the ef- 
fects of pairing and condensation by introducing non- 
condensed pairs appropriate for high effective magnetic 
fields. First, we introduce a Landau level representa- 
tion for the Gor'kov equations. The real-space Gor'kov 
coupled equations for the gap A(r) and the fermionic 
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Green's function G(r,r f ;iu) are: 

x A(r")G°(r'",r";-iuj)A\r'")G(r'",r';iuj) (1) 
A f (r) = 7jJ2f dr'G{r',r;iLu)G°{r',r;-iuj)A\r') 

iu; 

(2) 

Throughout this paper iu (id) will be used to de- 
note discrete fermionic (bosonic) Matsubara frequen- 
cies, with the traditional subscripts omitted for clar- 
ity. Introducing a Landau level basis for the fermions 
indexed by m = (N,p,k z ) where TV is the Lan- 
dau level, p the degenerate Landau level index, and 
k z the momentum parallel to the magnetic field, 
we write the bare Green's function G°(r,r f ;iu) = 
Z) n ^n(r)^ n {r r ) / {iuj - i n ) where £ n is the single-particle 
energy; the dressed Green's function is G(r,r';iu) = 
E^^'M^W^M [HI. The self energy, 
given by T,(r,r']iuj) = — A(r)A^ (r')Go(r' , r; — iu), is 
rewritten as Z mrn ,(iu) = - ]T n G n (-icj)A mn A] n/n , and 
the number equation necessary for a self-consistent so- 
lution is TV = ^Em^^mmM. Defining A mn 
J <irA(r)'0j n (r)'0^(r), and integrating over position vari- 
ables yields 

G mmj >{iu) = G^(zu;)<J mm / 

-^G:HA ml Gf(-^)A} n G nm 'H- (3) 

In 

'iEE^-(-«H4 (4) 

iuj mm'n J I \ / 1 

These correspond to Eqs. Q-Q respectively in the Lan- 
dau level basis. E q. Q is equivalent to that found else- 
where in Refs. [TSJ I16| 7 As in these references we restrict 
our consideration to intra-Landau level pairing. This as- 
sumption is important for arriving at a tractable scheme 
and a good approximation in the high-field regime [17] . 

It is useful to recognize that the nonlinear gap equa- 
tion, Eq. Q, applies to all T < T c . This nonlinearity 
is reflected in the presence of one dressed and one bare 
Green's function, as opposed to the two bare Green's 
functions associated with the instability onset in strict 
BCS theory. This gap equation can be thought of as a 
Bose-Einstein condensation condition [18 which reflects 
the vanishing chemical potential of the pairs below T c , 
and which becomes finite above. 

Characterizing Non-condensed Pairs Outside of 
the weak-coupling limit, pairs may form in kinetically 
excited states [18]. These pairs appear above the su- 
perfluid transition temperature T c , corresponding to a 
pseudogap phase, and should also persist below the tran- 
sition. To include these non-condensed pairs, we must 



therefore introduce variables corresponding to their exci- 
tation parameters. 

In the z-direction parallel to the field, we introduce a 
total momentum q z = k Zl + k Z2 of the pair, in analogy 
with the zero-field case [18]. Perpendicular to the effec- 
tive magnetic field, however, both condensed and non- 
condensed pairs may lie in the same Landau level. We 
posit, therefore, that non-condensed pairs correspond to 
different real-space gap configurations A(r), which we 
now parametrize by f. We single out Co which denotes the 
condensate gap configuration. In general non-condensed 
pairs can occupy other functional forms of A(r), for 
which the parameter ( will generally be close to £o, asso- 
ciated with low-energy excitations as described below. 

In our mean-field approach, we consider the conden- 
sate function A(r, Co) to be that of the optimal triangu- 
lar lattice. Non-condensed pairs can then occupy a two- 
dimensional continuum of other Abrikosov lattice con- 
figurations. Mathematically, we use the Landau gauge 
A = (0,Hx,0) and an Abrikosov lattice with unit vec- 
tors a = (0,a,0) and b = (b x ,b y ,0) : where ab x = 
with Ih = \/hc/eH the magnetic Hall length. The two 
clear mean-field distortions available to excited pairs are 
associated with changing b x /a, and changing b y /a (see 
Fig. [I]) [19 . Both of these are higher in energy rel- 
ative to the optimal Abrikosov lattice. We then fol- 
low Ref. [19] in setting ( = b y /a + ib x /a, for which 
the optimal triangular lattice condensate configuration is 
Co = l/2+iv // 3/2. To transform to state-space, we can as- 
sociate with each distortion a normalized real-space gap 
configuration A°(r,C), where J dr |A°(r,C)| = 1, from 

which A^ n (C) = J drA°(rX)^ti( r ) rl l J U r ) can also be 
calculated. 

In summary, these non-condensed pairs are now spec- 
ified by three degrees of freedom, thereby compensating 
the problematic dimensional reduction and restoring the 
possibility of stable condensation. We can also use these 
parameters to rewrite the gap equation, Eq. Q, by defin- 
ing a pair susceptibility x(C Qz] ify = 

-a J2 ^mm'(0 G rn^(i^)G%(q z -k z ;iQ-iu) (5) 

so that Eq. Q assumes a simple and suggestive form 

l+3X(Co,0;0) = 0. (6) 

In Eq. ^, G° is written in terms of the Landau level 
N and z- momentum of m = (iV,p, k z ) (with N m = N m > 
and k Zm = k Zml ) and <f? mm ,(Q = £„ A^ n (C)A°L'(C)- 
Green's Function of Non-condensed Pairs 
Equation (|6| suggests that there is a t-matrix (or a 
summation of particle-particle ladder diagrams), which 
is related to the pair susceptibility Eq. ([5| and which 
diverges at and below T c . More precisely, we have shown 
in Ref. [20] that, not surprisingly, this t-matrix is given 

by 

*p«(C, q z ; in) = — — 9 - — . (7) 
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FIG. 1. (Color online) A diagram of the different distortions 
of the Abrikosov lattice, the resulting splitting in energy lev- 
els, and the different occupation statistics that result, (a): 
The shaded gray unit cell is the optimal lattice configuration, 
whereas the red dotted unit cell corresponds to an excitation 
of by /a while the blue dashed unit cell corresponds to an ex- 
citation of b x /a. (b): Real-space diagrams of three different 
values for b x /a, showing zeroes of A°(r) for each configura- 
tion (black circles) and lattice vectors (dashed red lines). Eo 
is the optimal configuration, while E\ and E2 are progres- 
sively higher in energy, (c): The pairs (here pairs of blue 
(spin up) and red (spin down) fermions) are now able to oc- 
cupy a continuum of energy levels corresponding to different 
lattice configurations. For the three displayed configurations, 
(i): The BCS approach results in only the optimal configura- 
tion being occupied (also the case in this system at T = 0). 
In contrast, in this system for higher energy levels can 

also be occupied by pairs, (ii): An example of occupation 
statistics at a low temperature - most pairs are in Eo. (iii): 
At a higher temperature, more pairs are in excited states. 



Moreover, this t- matrix will lead to a self energy con- 
tribution [20 in parallel with what is found in Gor'kov 
theory via the condensate t-matrix, X^ m/ (ioj) = 



4> 2 mm >(Ot pg ((, q z \ in)G° N (q z -k z ;itl- iw). (8) 



This self energy is then fed back into the Green's function 
which enters into the pair susceptibility, and into the self- 
consistently determined t-matrix. 

Because of the mixing through the self-energy of an 
infinite number of real-space gap configurations, the cal- 
culation of the pair susceptibility x 1S n °t analytically 
tractable. To make progress without the distraction of 
heavy numerics, we approximate the pair susceptibility 



as ~ 



f (C;iuj)G° N (q z -k z ;in-iuj) (9) 



too m,m 



where we have, in effect, decomposed G mrn > into sepa- 
rate contributions each associated with a distinct lattice 
structure. In this way Eq. ([8]), which becomes 



,((;iuj) « -^(() 



x G° N (q z 



iQ — iuj). 



(10) 



determines the Green's functions G mrn '(£;iuj). 

Importantly, we can further write that Jj^ m ,((]ioj) « 
- 0mm'(O |A pg | 2 G° N (-k z ; -iu), where the gap contribu- 
tion from the non-condensed pairs is given by |A pg | 2 = 
— q z ifi^(Cj^j^)- To derive this simplified (BCS- 
like) expression for the pseudogap self energy in a strong 
"magnetic field" we have used the fact that the small 
chemical potential of the pairs implies that £ pg is strongly 
peaked around (q z ,iQ) — (0,0) near and below T c . 
When we consider a diagonal pairing scheme [16], the 
Green's function that results has a familiar BCS form 
G mm ((,iuj) + 



M 2 -e 2 -iA^re m (c)_. 

In effect, the main approximation we have made is to 
consider each vortex configuration as an "independent 
system," sharing self-consistently the same magnitude 
of the energy gap |A pg | 2 , but possessing a distinct real- 
space gap parameter A°(r,£) which in turn provides a 
unique form factor ^> 2 nm /(C) i n the self-energy. This is 
similar in approach to considering an LG energy func- 
tional for the different forms of the gap parameter. In 
contrast to LG theory, here we explicitly incorporate the 
fermionic nature of the pairing. 

Fermionic Constituents of the Pairs This discus- 
sion has been cast in terms of the basis eigenstates ip m (r) 
which we now need to specify. Here we use the mag- 
netic translation group basis p~6j [2TJ [22] , although we 
have also explored an alternative orbit-center based pair- 
ing basis p~5j [20]. Magnetic translation group (MTG) 
pairing relates to the Abrikosov lattice and is associated 
with a Bloch-like index k = (k x ,k y ). The unit cell for 
the MTG has unit vectors 2a and 6, which means in 
turn that the state basis is dependent on £. Pairing 
for the £ appropriate to the basis occurs with a sin- 
gle partner between opposite k [16], i.e. qz ( r ) = 

^N™L° n t(^)^N™i C ;- fc2+9 ,,i(^)- A^„(C) for these pairs is 
computed elsewhere [I6j [20] . 

Results With this formalism in place, we are in a po- 
sition to examine the underlying physics of how super- 
fluidity with pre-formed pairs takes place in the presence 
of a high effective magnetic field. This addresses the 
difficulty that a stable transition seems to require some 
inhomogeneity in the normal state [8] . In Fig. [5] we plot 
the t-matrix vs. lattice configuration £ for three differ- 
ent sets of effective temperatures, demonstrating that as 
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FIG. 2. (Color online) A three-dimensional plot of |£ pg (C,0;0)| against lattice configurations f. These plots are slightly above 
T c , such that in the rightmost plot — x(0>0;0) — — Mpair/^o where /x pa ir/^o is chosen to be —3 x 10 -4 (see Ref. [18] for a 

description of these parameters), and the other two plots share the same interaction g and fermionic chemical potential fi — 1.0, 
but different (A pg ,T): (a) (5.75,1.1118); (b) (5.90,1.0480); (c) (6.00,1.00). The fermions lie in the lowest Landau level for 
simplicity. 
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FIG. 3. (Color online) (a): A density plot of the total 
squared energy gap |A PS | 2 in real space, corresponding to 
Mpair /Z = -3 x 10~ 4 , fj, = 1.0, A pg = 6.0, and T = 1.0 
(the same parameters as in Fig. [2J;). This calculation is done 
through a discrete sampling of 1,360 points in C-space, and 
normalized to the largest value of the gap. (b): A density 
plot of the condensate energy gap |A(r, Co)| 2 only, which cor- 
responds to the total energy gap at zero temperature, again 
normalized to its largest value. 



condensation is approached from higher temperatures, 
the occupation of lattice states near the ideal triangu- 
lar Abrikosov lattice (£ = (o) begins to peak. Precisely 
at T = T c , a delta function results at C = Co- Never- 
theless at the transition there is still considerable weight 
associated with other lattice configurations reflecting the 
fact that the pseudogap |A pg | 2 remains finite. The con- 
densate contribution corresponds to a perfect triangular 
lattice which necessarily has small weight near T c . 

Of particular interest is the real space reflection of 
these distorted Abrikosov lattice contributions. To il- 
lustrate this precursor vortex configuration we evaluate 

l APS WI 2 = ^ E ^ ps (C,^;^)|A°(r,c)| 2 , (11) 



which is a weighted average of the gap (squared). This 
is compared with the counterpart for a fully condensed 



system in Fig. [3] By addressing the square of the gap, we 
emphasize that there is no phase information in the nor- 
mal state pseudogap. It should be noted that the point 
r = is chosen as a point of "symmetry breaking" or 
pinning center which breaks the translational symmetry 
available in the selection of each A°(r, Q. 

Conclusions In this paper we have addressed an im- 
portant and in principle testable prediction in the cold 
Fermi gases: the presence of precursor vortex configura- 
tions in the normal (pseudogap) phase. Indeed, the role 
of non-condensed bosons in rapidly rotating condensates 
has not been elucidated even for the atomic Bose gases. 
Here, too, one might expect a precursor vortex configura- 
tion. We have, moreover, elucidated the nature and role 
of excited pair states throughout the BCS-BEC crossover 
in high effective magnetic fields, showing that these are 
associated with distortions of the Abrikosov lattice. 

The concept of a "normal state vortex" liquid is also 
rather widely discussed in the context of high T c super- 
conductors [2[ [3]. There are, of course, a host of con- 
troversial issues associated with the cuprate pseudogap, 
but Fig. [3] illustrates one scenario for how one might think 
about this phenomenon in the context of pre-formed pairs 
in a strong magnetic field. 

The experimental implementation of this work in cold 
gases will depend on reaching a regime in which the 
Landau level spacing is much larger than the gap, al- 
though we expect the qualitative aspects regarding the 
role of pair density inhomogeneities to apply to much 
lower fields. While rapid rotation may be able to reach 
this regime, current proposals for artificial fields [1H6] also 
show promise. In this paper, two spin states which in- 
teract via a Feshbach resonance must both feel the same 
effective field, which may make geometric gauge-based 
proposals more challenging. Regardless of the implemen- 
tation, reaching this regime would facilitate a large body 
of new research into BCS-BEC phenomena. 

We are grateful to Jonathan Simon for helpful discus- 
sions. This work is supported by NSF-MRSEC Grant 
0820054. P.S. acknowledges support from the Hertz 
Foundation. 
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